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ABSTRACT 

Aims. To study cosmological models in scalar tensor theories of gravity with power law potentials as models of an accelerating universe. 
Methods. We consider cosmological models in scalar tensor theories of gravity that describe an accelerating universe, and we study a family 
of inverse power law potentials, for which exact solutions of the Einstein equations are known. We also compare theoretical predictions of our 
models with observations. For this we use the following data: the publicly available catalogs of type la supemovae and high redshift Gamma 
Ray Bursts, the parameters of large scale structure determined by the 2-degree Field Galaxy Redshift Survey (2dFGRS), and measurements of 
cosmological distances based on the Sunyaev-Zel'dovich effect, among others. 

Results. We present a class of cosmological models which describe evolution of homogeneous and isotropic universe filled in with dust like 
matter and a scalar field that is non minimally coupled to gravity. We show that this class of models depends on three parameters: Vo - the 
amplitude of the scalar field potential, Hg - the present value of the Hubble constant, and a real parameter s which determines the overall 
I ■ evolution of the universe. It turns out that these models have a very interesting feature of producing in a natural way an epoch of accelerated 
J — ' expansion. We fix the values of these parameters by comparing predictions of our model with observational data. It turns out that our model is 
compatible with the presently available observational data. 

Key words, cosmology: theory - cosmology: quintessence - large-scale structure of Universe — Noether symmetries-Scalar tensor theories 

1. Introduction 

Recent observations of the type la supernovae and CMB anisotropy indicate that the total matter-energy density 
of the universe is now dominated by some kind of dark energy causing an accelerated expansion of the Universe 
(Perlmutter 1997, Riess & al. 1998, Riess & al. 2004, Spergel & al. 2006). The origin and nature of this dark energy re- 
mains unknown (Zeldovich 1967, Weinberg 1989). 

Prompted by this discovery recently a new class of cosmological models has been proposed. In these models the standard cos- 
mological constant A-term is replaced by a dynamical, time-dependent component - quintessence or dark energy - that is added 
to baryons, cold dark matter (CDM), photons and neutrinos. The equation of state of the dark energy is assumed to be of a hydro- 
dynamical type W0 = P4,/p0, where and are, respectively, the energy density and pressure, and -1 < < 0, what implies 
a negative contribution to the total pressure of the cosmic fluid. When = -1, we recover the standard cosmological constant 
term. One of the possible physical realizations of quintessence is a cosmic scalar field (Caldwell, Dave & Steinhardt 1998), 
which induces dynamically a repulsive gravitational force, that is responsible for the observed now accelerated expansion of the 
universe. 

The existence of dark energy that now dominates the overall energy density in the universe is posing several theoretical problems. 
Firstly, it is natural to ask, why do we observe the universe at exactly the time when the dark energy dominates over the matter 
(cosmic coincidence problem). The second issue, a fine tuning problem, arises from the fact that if the dark energy is constant, 
like in the standard cosmological constant scenario, then at the beginning of the radiation era its energy density should have 
been vanishingly small in comparison with the radiation and matter component. This poses a problem, since in order to explain 
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the inflationary behaviour of the early universe and the late time dark energy dominated regime, the dark energy should evolve 

and cannot simply be a constant. All these circumstances stimulated a renewed interest in the generalized gravity theories, and 
prompted consideration of a variable A term in more general classes of theories, such as the scalar tensor theories of gravity. 

In our earlier paper (Demianski & al. 2006) we have analyzed extended quintessence models, for which exact solutions 
of the Einstein equations are known, and discussed how in such models it is possible to treat the fine tuning problem in an 
alternative way. We applied our consideration to a special model, based on one of the most commonly used quintessence potentials 
V(<p) - A<p'^, corresponding to the coupling Fief)) - (3/32)0^ (so called induced gravity). We showed that this model corresponds 
to a special, and physically significant case, which emerged by requiring the existence of a Noether symmetry in the pointlike 
Lagrangian. In this paper we analyze a new and wider class of theories derived from the Noether symmetry requirement. One 
of the main advantages of such models is the fact that they exhibit power law couplings and potentials, and admit a tracker 
behaviour. In some sense we complete and generalize the analysis initiated in (Marino & de Ritis 2001, de Ritis & al. 2000), 
where the attention was focused on the mechanism of obtaining an effective cosmological constant through the cosmological no 
hair theorem, and the analysis of the solution was restricted to the asymptotical f ^ oo regime only. Extending our analysis to 
the whole time evolution, we are able not only to clarify the properties of such solutions, but also to compare predictions of such 
models with observations. We concentrate on the following data: the pubUcly available data on type la supemovae and Gamma 
Ray Bursts, the parameters of large scale structure determined by the 2-degree Field Galaxy Redshift Survey (2dFGRS), and the 
measurements of cosmological distance with the Sunyaev-Zel'dovich eff'ect. 



2. Model description 

2.1. Specifying ttie model 

Since the detailed properties of a quintessence model, whose coupling and potential form is derived by requiring the existence 
of a Noether symmetry, are extensively discussed in (Demianski & al. 2006, from this time on Paper I), here we only summarize 
the basic results, referring readers to our previous paper for details. Let us consider the general action for a scalar field (f>, non 
minimally coupled with gravity, but not coupled with matter, in this case we have 

^ = ^ [Pim + ^^''V.^.^.v - V(0) + Xmj d^X, (1) 

where F {</>), Vi(f>) are two generic functions representing the coupling of the scalar field with geometry and its potential energy 
density respectively, R is the scalar curvature, -g^^<f>,jj^,v is the kinetic energy of the scalar field (f>, and £,„ describes the standard 
matter content. In units such that SnGff -h-c- 1, where Gn is the Newtonian constant, we recover the standard gravity when 
F is equal to -1/2, while in general the effective gravitational coupling is G^// = Here we would like to study the simple 
case of a homogeneous and isotropic universe, what implies that the scalar field depends only on time. It turns out that for the 
flat Friedman-Robertson- Walker universe filled with matter satisfying the equation of state p = {j — V)p and the scalar field 0, 
the action in Eq. (1) reduces to the pointlike Lagrangian 

£, = 6Faa^ + 6F'4>a^a + c? - V(0)j - Da"^*^"^' , (2) 

where a is the scale factor and prime denotes derivative with respect to 0, while dot denotes derivative with respect to time. 
Moreover, the constant D > is defined in such a way that the matter density is expressed as pm = D{aoldf'^, where 
1 < 7 < 2. The efl'ective pressure and energy density of the 0-field are given by 

P0 = ^<^'-V(0)-2(F + 2//F), (3) 

P^ = \ip'' + V{4>) + 6HF. (4) 

These two expressions, even if not pertaining to a conserved energy-momentum tensor, do define an effective equation of state 
p,p = w^p^, which drives the late time behavior of the model. 

From now on we restrict ourselves to a dust filled universe with y - I, and - 0. Using the point like Lagrangian Eq. (2) 
in the action and varying it with respect to (p we obtain the Euler-Lagrange equations 

ip + 3Hip + 6{H + 2H^)F' + V' = 0, (5) 
2H + 3H^=^{p^ + Pm), (6) 
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together with the first integral 

H'^ = -^{p^+p„) . (7) 

Let us now introduce the concept of an effective cosmological constant h^fj. Using Eq.(7) it is natural to define the effective 
cosmological constant as K^ff = and the effective gravitational constant as G^ff - -jp- With this definitions we can rewrite 
Eq.(7) as 

3H^ = GeffPm + Kff. (8) 

Introducing the standard Omega parameters by 
o P« o _ ^«// _ P'l' 

we get as usual that 

a„ + Qa,^^ = 1 . (9) 

Imposing the Noether symmetry in the quintessence minisuperspace, {a, (p; d, 0}, where the point-like Lagrangian is defined, it is 
possible to exactly integrate the Einstein field equations Eqs. (7) and (6), as well as to find a form for the two unknown functions 
F((p) and V((p)) (for details see, Capozziello & al. 1996, Demianski & al. 1991). The existence of this symmetry actually leads to 
the following relation between the functions F{(/>) and V((py. 

V ^ VQ(F(4>)r^'\ (10) 

where Vq is a constant and 
3(s+ 1) 

with s a real number. Moreover, a possible simple choice of the coupling is 

F = mi<f> + <Po)\ (12) 
where (f>o is a constants that does not affect our results and, therefore, from now on we set it to zero, and 

^^'^ = 48(.-H)(. + 2) • ^''^ 

Let us note that the form of the coupling given by (12) is quite relevant from the point of view of fundamental physics. It 
describes the so called Induced Gravity. The Induced Gravity model was initially proposed by Zee in 1979 (Zee 1979), as a theory 
of gravity incorporating the concept of spontaneous symmetry breaking. It was based on the observation that in gauge theories 
the dimensional coupling constants arising in a low-energy effective theory can be expressed in terms of vacuum expectation 
values of scalar fields. In such a model the gravitational and cosmological constants are not introduced by hand, but are generated 
in the process of symmetry breaking of a scalar field non minimally coupled with the Ricci scalar in the Lagrangian describing 
the system. Once the Noether symmetry is specified it is possible to find a corresponding conserved quantity and use it as a 
new dynamical variable (for details see (Capozziello & al. 1996, Demianski & al. 1991)). One can then solve the corresponding 
Lagrange equations and finally after returning to the original variables we obtain the sought after a{t) and (p{t). The final results 
can be written in the form 

a(t) = A(s) lB(s)t~ + —1 , (14) 

2.T+3 

<pit)^Cis){-^Bis)t- + —] t-—>, (15) 
where A(i), B{s), C{s), y{s) and;^(i) are given by 

Ms) - ixis))-[-^] , (16) 
Bis) = I^^M^^Y"^ , (17) 



s + 6 / \ 3y{s) 

s + 3): 



_(3+2£) 

(2«3) (s + 3)2o\ 

Cis) = (xis))-—\ \ // , (18) 
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Fig. 1. Diagram of ^(s) as function of s. It turns out tliat an attractive gravity requires s e (-2 , -1). 



and 

25 + 3 
12(s + i){s + 2) 

2i + 3 

where D is the matter density constant, Eq is a constant of integration resulting from the Noether symmetry, and Vo is the constant 
that determines the scale of the potential. Together with the independent parameters, we then use these three constants {D, So, Vo), 
which however are not directly measurable, but they parametrize the possible solutions of the model. In the next section we shall 
reduce their number by means of additional assumptions along the lines of Paper I. 



2.2. The parameter space 

As it is apparent from Eq.(14) and Eq.(15) there are two additional particular values of s, namely s - and i = -3 which should 
be treated independently. 

When i = 0, a Noether symmetry exists if: 

1. F ^ pQcfP-, and V = Vo^^ , 

2. P - ~2 (i^iiiii^^l coupling), and V - VoiAe^"^ - Be''^'^)^, with fj. = V3/2, and A, B being constants. 

The case of the minimal coupling has been thoroughly investigated in (Rubano and Scudellaro 2002, Rubano & al 2003, 
Demianski & al. 2005). If B = we obtain an exponential potential, which is a very important model of quintessence with a 
standard scalar field. When i = -3 we recover the case of the quartic potential treated in the Paper I. We shall therefore con- 
centrate on the other values of s. As we shall see in a moment, this will lead to a very different class of potentials from those 
discussed in Paper I. In fact we obtain inverse power-law type potentials, which are interesting and recently widely used in the 
literature. 

First of all, we have to find the physically acceptable range for s, and the most important requirement is, of course, that 
Geff > 0, i.e. F <Q. This restricts s to (-2,-1) as shown in Fig.(l). 

As mentioned above in the range s € (-§ , -1) the potential V((/)) is of an inverse power-law type, In this case 

our model naturally admits cosmological scaling solutions, recently studied, for example by (Amendola 1999, Uzan 1999) in the 
context of quintessence models. In Fig. (2) we see that all the possible exponents for the inverse power-law potential are available. 
To determine the integration constants D, and So we follow the procedure used in Paper I, and we use the age of the universe, 
fo, as a unit of time. Because of our choice of time unit, the expansion rate H{t) is dimensionless, so that our Hubble constant is 
clearly of order 1 and is not (numerically) the same as the Hq that is usually measured in kms^'Mpc"'. Setting oq - a(to) = 1 
and //o - //(fo), we are able to express So and D as functions of s and Hq. We obtain: 



D 



^ ■ B(s) 



,Ais) 



(.^+l)(^+3) 

3.t-+7.T+3 / _ . ^ _ _ _\ \ s^-s-i 



3-^^^(s + 3y^^^(s + 6) [(Ho - 2)s^ + 3(//o - 2)s - 3)y(s)7TT. 



(21) 



(22) 




-20 



Fig. 2. Behaviour of the couphng factor ^(i) (red curve) and the power-law exponent p(s) (blue curve): we see that with an appropriate choice 
of s in the range (-1.5 , -1) all the values for the exponents are available. 
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Fig. 3. Plot of log[oP^ versus log,Qa (solid black line). The upper and lower dashed lines indicate the log-log plot of a" 
respectively. It turns out that scales as a"", with 3 < n < 4. In this and subsequent plots we use the mean values for the parameters obtained 
through fits (see Table 2). 

Therefore, our family of models, labeled by s, depends only on Vo, and //q. For both these parameters we have not only a thorough 
knowledge of their physical meaning, but for Hq we can also strongly constrain its range of variability. Actually, we may easily 
obtain the relation 

h = 9.9—, (23) 

T 

where as usual h = //q/IOO and r is the age of the universe in Gy. We see that Hq fixes only the product hr. In particular, we 
know that t - 13.73^q [j (see for instance Spergel & al. 2006), thus for //q « 1 we get h < 0.76. The actual value of h may be 
obtained by some of the subsequent tests or in others it has to be set as a prior - this will be specified in each case considered 
below. 

Using the available observational data we can further constrain the range of possible values of s. Actually requiring that 
today a(tix) > 0, as indicated by observations of supernovae la and WMAP, we constrain the range of possible values of s to 
i e (-1.5, -1.2). 

From the Eqs. (14) and (15) it turns out that for small values of t the scale factor and the scalar field behave as 

a oc t ■(■+3> , (24) 

(2.!+3)- 

(p OC t 2.'(.'+3) . (25) 



Substituting these functions in p^, as given by the Eq.(4), we get that for small f the scalar field density oc a^^. This is however 
true only asymptotically for very small t. Exact computation as shown in Fig. (3) gives oc a ", with n ~ 3. This justifies our 
assertion that our model naturally admits scaling solutions. The situation changes dramatically near the present time (see below). 
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Fig. 4. as a function oi x = Logioil + z), for the averaged mean values provided by our analysis, as 
shown in Table 2: we observe a transition from a small constant value in the past, Iw^l « 0, to vw^ = -1 
at present. 



For large values of t the scale factor and the scalar field behave as 



(26) 
(27) 




respectively, as shown in (Marino & de Ritis 2001), where this asymptotical regime is discussed. It is interesting to note, as 
shown in Paper I, that is representing an equation of state, in the usual sense, of the effective cosmological constant A^//, and 
that Aeff asymptotically behaves as a true cosmological constant (w^ -1) as / — > oo. Let us also mention that since both 
and depend on F((p) through its time derivative (see Eqs. (4) and (3)) and asymptotically F(0) ~ constant we recover in this 
limit the minimally coupled theory. In Figs. (4) and (5) we show the evolution with the redshift and the rate of evolution of wf, 
the fast transition toward — > - 1 appears between z = 3 and z - 0.5. Before reaching this asymptotic regime the energy density 

is dominated by the couphng term 6HF. Concluding this section we present the traditional plot logp^ - log a and compare 
it with the evolution of matter density (see Fig. (6)). It is interesting to see that follows the matter density during the matter 
dominated era, and it becomes dominant at late time. 

3. Newtonian limit and Parametrized Post Newtonian (PPN) beiiaviour 

Recently the cosmological relevance of extended gravity theories, as scalar tensor or higher order theories, has been widely 
explored. However, in the weak field approximation, all these classes of theories are expected to reproduce the Einstein general 
relativity which, in any case, is experimentally tested only in this limit. This fact is a matter of debate since several relativistic 
theories do not reproduce Einstein results at the Newtonian approximation but, in some sense, generalize them, giving rise, for 
example, to Yukawa-like corrections to the Newtonian potential which could have interesting physical consequences. Moreover, 
in general, any relativistic theory of gravitation can yield corrections to the Newton potential (see for example. Will 1993) which 
in the post-Newtonian (PPN) formalism could furnish tests for such theory, mainly based on the Solar System experiments. In 
this section we want to discuss the Newtonian limit of our class of scalar-tensor theories of gravity, the induced gravity theories, 
and to study the Parametrized Post Newtonian (PPN) behaviour of these theories. In particular, it turns out that the Newtonian 
limit depends on ^(s). Furthermore, we find a quadratic correction to the Newtonian potential strictly depending on the presence 
of the scalar-field potential which acts as a cosmological constant. 

3.1. Newtonian limit 

In order to recover the Newtonian hmit of our theory described by the action in Eq. (1), associated with the effective stress-energy 
tensor 




(28) 




Fig. 5. Rate of change of the equation of state as measured by versus the parameter. The values of the parameters correspond to the 
average values provided by our analysis and shown in Table 2. 



where T^'J, is the usual stress-energy tensor of matter and □ is the d' Alembert operator, we write the metric tensor as 

guv - rjfiv + hfiy , (29) 

where //^v is the Minkowski metric and /i^v is a small correction. In the same way, we define the scalar field as a perturbation, 
of the same order as h^y, of the original field (f), that is 

= ^0 + -A, (30) 

where ip^ is a constant of order unity. It is clear that for = 1 and i/' = the standard Einstein general relativity with A is 
recovered. 

To write the Einstein tensor G^y in an appropriate form, we define the auxiUary fields 

h^v = hf,y - ^T]^yh, (31) 

and 

O-a = haf},ylf^ , (32) 

where h = rf^h^y. Given these definitions, to the first order in h^y, we obtain 

Gf,v - [a,,h^y + rj^yCTafirf^ - cr^^y - o-y,^} , (33) 

where u^(p = //"'^jjy. To obtain the weak-field limit of the field equations we have also to expand the effective stress-energy 
tensor: this means that it is necessary to expand the coupUng function and the self interacting potential. Specifically, it turns out 
that expanding the coupling function F{(f>) and the self-interacting potential V(<p) (by using their explicit forms) up to the second 
order in tfr, we get 

Ficf>)^^(s){ipl+2ipoif,+ ifr") , (34) 

V(4>) ^ Vq Up'^; + p(s) (f'f^' ifr + (p^q' + ■ ■ ■ I , (35) 

where Vo - VQ^(sy^''\ Then, to the first order, the effective stress-energy tensor becomes 

T-, 2+p(s) 2 

V = -2<Po %ya„ilr + 2 (fiQi^^f^y 2ii7r'^'" ~ 2 ^ ^ ^ 

and the field equations assume the form 

1 [aji^y + ri^ya-afirj"'^ - cr^^y - o-^,,^} = 2 iplq^yU^if/ - 2 (filiff,f,y+ (37) 

^ 2^is) '^"^'^2 05)^'^'" 
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Fig. 6. Plot of logidP^ versus logjoO in the Jordan frame . The vertical bar marks logmflo- The solid red 
straight line indicates the log-log plot of p„, versus a. The matter dominated era and the transition to the 
present dark energy dominated regime are represented. 



When describes a point particle of mass M and for p(s) + 4 and p{s) ?i 1 we get 



■200 



1 - 16^(.) 



M 
r 



M^l %{p{s)-A){p{s)-l) 



1 - \2^{s) 



%{p{s)-A){p{s)-\) 



^(s) 2{p{s) - 1) 



hij ^ 6ij- 



4jt 



2 1 - 8^(^) 



M 

h 

r 



I -2m 

Mifl Voipis) - 4-)(p(s) - 1) 



1 - 12^^) 



l-2m 



f> 2+p(s) 



H Voipis) -4)(p(s)-l) 



m 2ipis)-l) I -2m 

where only terms linear in Vo are given and we omitted the constant terms. 



(38) 
(39) 
(40) 



3.2. Constraints on PPN parameters 

A satisfactory description of the PPN limit for scalar tensor theories has been developed in (Esposito-Farese 2004, 
Damour & al. 1993). In these papers, this limit has been thoroughly discussed leading to interesting results even in the case 
of strong gravitational sources like pulsars and neutron stars where the deviations from General Relativity are considered in 
a non-perturbative regime (Damour & al. 1993). The starting point of such an analysis is a redefinition of the non minimally 
coupled Lagrangian action in terms of a minimally coupled scalar field model via a conformal transformation of the form 
gfiv = - 2Fi(f>)g^y. In fact, assuming the transformation rules: 



(dcf>\^ _ 3 ld\nFi(P)\^ 1 



[d<p) 2\ d(p j 2Fi(/)) 
and 



(41) 



= = ^".^ (<^)' (42) 

4Fi(f>y 

one rewrites the action as 



R + ^?'>.;-0.v - Vi<f>) + £,r, 



d\. (43) 



The first consequence of such a transformation is that the non-minimal coupling is transferred to the ordinary matter sector, intro- 
ducing an interaction term between matter and the scalar field. Actually, the Lagrangian X.m depends not only on the conformally 
transformed metric and the matter fields but it also depends on the coupling function. In the same way, the field equations 

can be recast in the Einstein frame. The energy-momentum tensor is defined as T'" - -^=^^ and is related to the Jordan 

expression as T™, - ^-2 Ficp) TJJJ,. Possible deviations from the standard General Relativity can be tested through the Solar 
System experiments (Will 1993) and binary pulsar observations which give an experimental estimate of the PPN parameters. 
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The generalization of these quantities to the scalar-tensor theories allows the PPN-parameters to be expressed in terms of the 
non-minimal coupling function and in our case using the Eqs. (12) and (13), we obtain : 

/™_i = e:^ 

^ F(0) + 2[F'(0)]2 1 + 8^(5)' 

^ F(0) + 3[F'(<*)]2 # • 

The above definitions imply that the PPN-parameters in general depend on the non-minimal coupling function F((p) and its 
derivatives. However in our model depends only on s while ^^^'^ = 1 . The PPN-parameters can be directly constrained 
by the observational data. Actually, Solar System experiments give accurate indications on the ranges of ^ , jSg™. We 



Table 1. A brief summary of recent constraints on the PPN-parameters. 



Mercury Penh. Shift (Shapiro 1993) 


127,';™ -/^r"'-ii< 3x10-^ 


Lunar Laser Rang. (Williams & al. 1996) 


4y8o™ - To™ - 3 = -(0-7 ± l)xl0-3 


Very Long Bas. Int. (Shapiro & al. 2004) 


Ir^™ - 1| = 4x10-'* 


Cassini spacecraft (Bertotti, less & Tortora 2003 ) 


r^™ - 1 = (2.1 ±2.3)x 10-5 



sunnmarize the experimental results in Table 1. These results have been used by (Schimd, Uzan & Riazuelo 2005) to set the 
following constrains : 



TN 



11 < 2x10" 



11 < 6x10" 



(46) 



It turns out that the limit for yS^™ in the Eq. (46) is naturally verified, for each value of s, while the constraint on ly^™ - 1| is 
satisfied only for s e (-1.5, -1.4), as shown in Fig. (7). 

For the sake of completeness, here we take into account even the shift that the scalar-tensor gravity induces on the theoret- 
ical predictions for the local value of the gravitational constant as coming from the Cavendish-hke experiments. This quantity 
represents the gravitational coupling measured when the Newton force arises between two masses : 



^Cav = ■ 

OT1-OT2 

In the case of scalar tensor gravity, the Cavendish coupling is related to F and F' and is given by : 



Gcav - --. 



1 



2F(0) 



1 -I- 



[2F'(0)]2 



4F(0) + 3[2F'(0)]2 



1 



1 + 16«*) 



= G 



1 + 16^(5) 



1 -I- 12^(i) 



1-1-12^(5) 

and in our models it depends only on s 

Finally, in Fig. (8) we plot the Brans-Dicke parameter ubd as a function of s: actually, for our model u)bd = - 
1 



(47) 



(48) 



F(^) 



4^(5) 



2F'(0)2 

It turns out that just for s e (-1.5, -1.4), cjbd satisfies the Umits coming both from the Solar System experiments. 



o)bd > 40000 (Will 1993), and current cosmological observations, including cosmic microwave anisotropy data and the galaxy 
power spectrum data, give (Obd > 120 (Acquaviva & al. 2005).^ 



4. Some considerations about conformal transformations and interacting darl( energy 

Conformal transformations are often used to convert the non minimally coupled scalar field models into the minimally coupled 
ones, to gain mathematical simplification. The Jordan frame, in which the scalar field is nonminimally coupled to the Ricci 
curvature, is mapped into the Einstein frame in which the transformed scalar field is minimally coupled but at a price of coupling 
matter to the scalar field. The two frames are not physically equivalent, and some care has to be taken in applying this technique 
(see for instance (Faraoni 2000) for a critical discussion about this point). In this section we study the efl'ect of conformal 

' We indicate with the subscript that the estimates are based on Solar System observations. 

^ It is worth noting that the legitimacy of the procedure of direct comparison of local and cosmological observations, in order to esti- 
mate variations of the physical constants, is a rather strong assumption, which deserves a proof or at least a justification. There is actually 
no reason a priori why local experiments should reveal variations occurring on cosmological scales, and in regions which are participating 
in the Hubble expansion. This interesting aspect of scalar tensor theories of gravity is discussed in (Clifton, Mota & Barrow 2005) and in 
(Shaw & Barrow 2005), where it is shown that such procedure is correct. 




-1.8 -1.6 -1.4 -1.2 -1 



Fig. 7. Current limits on tlie PPN parameters restrict the range of the parameter s. We see that the constraint on ly^™ - 1| < 2 • 10 leads to 
s € (-1.5,-1.4), as shown in the inner zoom. 
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Fig. 8. Behaviour of the Brans-Dicke parameter cobd as a function of s. For s e (-1.6,-1.3), cobd satisfies limits placed by the solar system 
experiments {cobd > 40000) and by current cosmological observations (cobd > 120) 



transformations on our models and show that, in presence of matter, it can mimic a coupling between the quintessence scalar 
field and dark matter. We discuss some implications of such a fictitious interaction on the effective equation of state. Actually 
it turns out that, since the interaction alters the redshift-dependence of the matter density, it is possible to obtain an effective 
transformed dark energy equation of state of < - 1 . Let us start from the transformation rules connected with the conformal 
transformation gy,v - -'^F{4>)giiv- 



d(f> V 2F(0)2 ' ^ ' 



a = ^j-2F(f)a, (50) 

dt I 

Jt - ^''^ 

~ V(0) 



with these new variables the Lagrangian in the Eq.(2) becomes 

Z = 3af + a^i^t- y(0) j + I,„, (53) 
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Fig. 9. Time evolution of the transformed scalar field (f). 



where now dot denotes derivative with respect to t, and 

I„($,a) = V-2F(0)X„ = yj-2Fi4,)Da-^^^-'\ (54) 

An observer in the Einstein frame would infer that the scalar field is coupled to the dark matter, and this interaction is represented 
by the term fi(f>i(p)) = -y/-2 F((p), as it can be seen more clearly from the field equations: 

^.3H^=-V,-'^-^ (55) 

3^' = ilt- W)) + m -jp- a-\ (56) 
\^ / /(0o) 

where p^^ is the actual matter density (treated here as dust), 0o the value of today, and j^-^ = D. It turns out that the density 

__, _ / 

of dark matter does not evolve as a ^, but scales dispmcc —. Also the Klein-Gordon equation (55) differs from the standard one 

because of the last term on the right hand side. As shown in (Das et al. 2005) interaction between dark matter and dark energy 
could result in an effective equation of state p^s = WeffPeff > mimicking the ACDM model. Actually w^s is defined by the matter 
continuity equation 

^=-3H{\+Weff)p,s, (57) 
at 

where 

Peff=T=r — = 1 +P,^- (58) 

It can be shown that Weff is simply related to 

WefF = ^ , (59) 
I - X 

where 

,._^(/gl_lV (60) 



\/(0o) , 

Since x = today, one has vv^*^! - w"^\ which is greater than or equal to -1. It turns out (Das et al. 2005), however, that / 
increases in time, so that x > and it is possible to have w^s < - 1 in the past. In such a way interaction between dark matter and 
dark energy could generate a superquintessence regime, provided that the observer treats the dark matter as non-interacting, and 
ascribes part of the dark matter density to the scalar field, as is shown in Eq. (58). On the other hand an interacting dark energy 
component could hide the effect of a non standard gravity, in the Einstein frame, provided that it is considered as the physical 
one. 
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Fig. 10. Evolution with the redshift of in the Einstein frame. 

4. 1. Conformally transformed coupled dark energy 

Let us consider our nonminimally coupled model characterized by the functions F{</)) = ^{s)<p^ and V{<f>) = Vo(^(0))^**^ respec- 
tively. According to the rules in Eqs. (49 - 52) we obtain the following relations between the transformed and original dynamical 
quantities: 




hi0, 



,s)~2 j.2p(s)-4 



(61) 

(62) 
(63) 



As we see from Eq. (62), in the transformed frame the cosmic evolution is mediated by the presence of the scalar field. The 
explicit form of Eq. (51), which coimects the cosmic time in both frames, can be written in analytical form, but is indeed rather 
complicated: actually it turns out that 



t = { 



2^^ +6^+9 

2C{s)s{s + 3)rT(^ 



1 \^ 



B(s)- 



VoB(s)t- 



,3 f 



-2s^-6s-9 
6s 



1+3 -2.^-6.-9 ^ 1 

2s ' 6.? 



yis) 



VoB(s)t- 



1 -^p::^: 



B(.s) 



(2^2 + 6s + 9)({^) - B(.)) i^-y(s) (3^) + VoBis)t^3 + B(i)r(^)) 



(64) 



From the Eqs. (61-64) we can evaluate the scalar field energy density and pressure, and the equation of state according to the 
usual definitions. 

As in the Jordan frame, it turns out that also in the Einstein frame vv^ ^ -1 for ? — » 00 (see for instance Fig. (10)). 

Concluding this section we present the traditional plot logp^ - log and compare it with logp„, - log a relation (see Fig. (11)). 
Interestingly we see that, just because of the interaction term, does not track anymore the matter during the matter dominated 
era, as it happens in the Jordan frame (see Fig. 6), but becomes dominant at earlier times. It is interesting to write down the 
efl'ective equation of state We/f (see Eq. 57), which mimics a CDM model. Actually in our case Wg/f is 



Weff = 



W0 

1^ 



Pino 

a^p^ \2^o - 1 



(65) 



It turns out that in our case Wg/f and are practically indistinguishable, since in our model the transformed scalar field (p weakly 
evolves with time as shown in Fig. (9) 



4.2. Non minimally coupled quintessence and mass varying neutrinos through conformal transformations 

In this section we briefly discuss how the scalar tensor theories of gravity could be involved in a cosmological model with mass 
varying neutrinos that mimic the dark energy, a quite diff'erent theoretical scenario of evolution of the universe that recently has 
been suggested by (Fardon & al. 2004). Let us recall that recently the mass differences between neutrino mass eigenstates (wi. 
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Fig. 11. Plot of log,oP0 versus logmT? in the Einstein frame. The vertical bar marks log,Qao- We see that, 
just because of the interaction term, does not track anymore the matter during the matter dominated 
era, as it happens in the Jordan frame (see Fig. 6), and becomes dominant at later times. 

m2, irii) have been measured in oscillation experiments (Lesgourgues and Pastor 2006). Observations of atmospheric neutrinos 
suggest a squared mass difference of Arn^ ~ 3 x IQ^^eV^, while solar neutrino observations, and results from the KamLAND 
neutrino experiment, point towards Am^ ~ 5 x lO^^eV^. While only weak constraints on the absolute mass scale {"Liny - 
mi +1712 + niT,) have been obtained from single j6-decay experiments, the double beta decay searches from the Heidelberg-Moscow 
experiment have reported a signal for a neutrino mass at > Act level ( Klapdor-Kleingrothaus and al. 2004), recently promoted to 
> 6cr level (Klapdor-Kleingrothaus 2006). This last result translates into a total neutrino mass ofLniy > l.leV at 95% c.l., but this 
claim is still considered as controversial (see Elliott and Engel 2004). It is known in the literature (Lesgourgues and Pastor 2006) 
that massive neutrinos can be extremely relevant for cosmology as they leave key signatures in several cosmological data sets. 
More specifically, massive neutrinos suppress the growth of fluctuations on scales below the horizon scale when they become 
non relativistic. Current cosmological data have been able to indirectly constrain the absolute neutrino mass to 'Eniy < 0.75 eV at 
95% c.l. (Spergel & al. 2006), and are challenging the Heidelberg-Moscow claim. However, as first noticed by (Hannestad 2005), 
there is some form of anticoiTelation between the equation of state parameter w and Swy. The cosmological bound on neutrino 
masses can therefore be relaxed by considering a dark energy component with a more negative value of w^ff than a cosmological 
constant. Actually it has been proved that the Heidelberg-Moscow result is compatible with the cosmological data only if the 
equation of state (with w being constant) is w^jf < - 1 at 95%. This result suggests an interesting link between neutrinos and dark 
energy (see for instance Brookfield & al. 2005 b, Brookfield & al. 2005 a, Bi & al. 2004, Kaplan & al. 2004, Amendola 2004). 
According to this scenario the late time accelerated expansion of the universe is driven by the coupling between the quintessential 
scalar field and neutrinos. Due to such coupling the mass of neutrinos becomes a function of this scalar field. Since the scalar 
field evolves with time the mass of neutrinos is not constant (mass-varying neutrinos): the main theoretical motivation for such 
connection relies on the fact that the energy scale of the dark energy is of the order of the neutrinos mass scale. Moreover, as 
discussed above, in interacting dark energy models, the net effect of the interaction is to change the apparent equation of state of 
the dark energy, allowing a so called superquintessence regime, with We/f < -1. Interestingly enough, if the Heidelberg-Moscow 
results are combined with the WMAP 3-years data and other independent cosmological observations, as the ones connected with 
the large scale structure - coming from galaxy redshift surveys and Lyman-a forests, or with the SNIa surveys, it is possible to 
constrain the equation of state to -1.67 < We/f < -1.05 at 95% c.l., ruling out a cosmological constant at more than 95% c.l. 
(see De La MacoiTa et al. 2007). 

In the following we will discuss the coupling between neutrinos and dark energy from the point of view of the conformal 
transformation, according to the arguments outlined in the previous section, i.e. we will show that the neutrino mass and scalar 
field coupling can be interpreted as an effect of conformal transformations from the Jordan to the Einstein frame, just as it 
happens for the coupling between the dark energy and the dark matter. For our purpose the neutrinos can be either Dirac or 
Majorana particles, the only necessary ingredient is that, according to (Pardon & al. 2004), the neutrino mass is a function of 
the scalar field. In the cosmological context, neutrinos should be treated as a gas (Brookfield & al. 2005 b) and described by 
the collisionless distribution function fix', p', r) in the phase space (where r is the conformal time) that satisfies the Boltzmann 
equation. When neutrinos are collisionless the distribution function / does not depend explicitly on time. We can solve then the 
Boltzmann equation and calculate the energy density stored in neutrinos (/o is the background neutrino distribution function): 




(66) 



with - 



+ my(<p) 'cP', a is the scale factor and q' - a p' is the comoving momentum. The pressure is 




(67) 
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Fig. 12. Redshift dependence of the second derivative of the scale factor. The transition from a decelerating to an accelerating expansion occurs 
close to z ~ 0.5, as predicted by recent observations of SNIa Zi = 0.46 ± 0.13 (Riess et al. 2007, Riess & al. 2004). 

From these equations we derive that 

Py + 3H(py + py) = — ^4>(pv - 3pv) (68) 

(note that here the dot denotes the derivative with respect to r). The Klein Gordon equation for the scalar field reads 

'$+2H^+a-'^Z - -a^^-^(py - 3py). (69) 
d(p dcf> 

We see that Eq. (69) is formally equivalent to Eq. (55), but now the term on the right hand side describes coupling of the scalar 
field to the neutrino mass. As shown above the interaction between neutrinos and dark energy could result in an effective equation 
of state Weff, defined in Eq.( 57), while is defined in the standard way by = w^^. As in the previous section the Wgff 
parameter is related to by 

Weff = (70) 
1 - X 



with 

- PvO 



my{<l)) _ ^ 



(71) 



Also in this case w^ff can be less than -1, as was pointed out in the context of models with dark matter/dark energy interaction. 
This circumstance could lead to an observational test to establish which of the frames, Jordan or Einstein, is the physical one, 
since the coupling between the quintessential scalar field is provided by the function F{(f>), which should drive not only the dark 
matter evolution (see Eq. (55)), but also the neutrinos mass variation and the evolution of the gravitational constant Geff = -jp- 

5. Observational data and predictions of our models 

5. 1. Constraints from recent SNIa observations 

In this section we present results of fits of predictions of our model to the best SNIa data sets presently available. As a starting 
point we consider the sample of 182 SNIa compiled in (Riess et al. 2007), which includes the 21 new Type la supernovae recently 
discovered with the Hubble Space Telescope (HST), and combines previous SNIa data sets, namely the Gold Sample compiled 
in (Riess & al. 2004) supplemented by the SNLS data set (Astier & al. 2005). Following the procedure described in Paper 1, we 
perform a analysis comparing the redshift dependence of the theoretical values to the observational estimates of the distance 
modulus, p - m - M, which in scalar tensor theories of gravity takes the form 

OT - M = 5 log Diiz) + 25 + AMg . (72) 
Here the presence of the correction term 

AMG=^log(^), (73) 
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Fig. 13. Observational data of the SNIa sample compiled by (Riess et al. 2007) fitted to our model. The solid curve is the best fit curve, compared 
with a standard ACDM model with Q.^ = 0.71 (red dashed line). 

describes the effect of time variation of the effective gravitational constant Geff on the luminosity of high redshift supernovae, 
and allows one to test the scalar tensor theories of gravity (Gaztanaga & al. 2001 , Uzan 2003) using the SNla data. Moreover, 
for a general flat and homogeneous cosmological model the luminosity distance can be expressed as an integral of the Hubble 
function as follows: 



where E(z) - is related to the Hubble function expressed in terms of z = ao/a(t) - 1. Let us note that the luminosity distance 
depends also on the Hubble distance ^ ( which does not depend on the choice of the unit of time). Snchfreedom allows us to fit 
h or the a priori unknown age of the universe t using the SNla dataset. We find that;^'^^^ = 0.98 for 182 data points, and the best 
fit values are Hq = 0.96+°{j^, s = - 1 .46+J}j;;^, which corresponds to ^a^^ = 0.69+[J{J^ We also get h = 0.72 ± 0.04. 
In Fig (13) we compare the best fit curve with the observational data sets. 

5.1 .1 . Dimensionless coordinate distance test 

After having explored the Hubble diagram of SNla, let us now follow the more general approach, suggested by Daly & Djorgovski 
(Daly & Djorgovski 2004), and already tested in Paper 1. Consider as a cosmological observable the dimensionless coordinate 
distance defined as : 



noting that y(z) does not depend explicitly on h, so that any choice for h does not alter the main result. Daly & Djorgovski have 
determined y(z) for the SNla in the Gold Sample of Riess et al. (Riess & al. 2004) which is the most homogeneous SNla sample 
available today. Since SNla allows to estimate Di rather than y, a value of h has to be set. Fitting the Hubble law to a large set of 
low redshift (z < 0.1) SNla, Daly & Djorgovski (Daly & Djorgovski 2004) have set: 

h = 0.66 + 0.08. 

To enlarge the sample, Daly & Djorgovski added 20 points on the y{z) diagram using a technique of distance determination 
based on the angular dimension of radiogalaxies (Daly & Djorgovski 2004). This data set has been recently supplemented by 71 
new supernovae from the Supernova Legacy Survey of (Astier & al. 2005), which allowed the determination of dimensionless 
coordinate distances to these supernovae. These were obtained using the values and uncertainties of yu^ listed in Table 9 of 
(Astier & al. 2005), with h = 0.7. This extended sample that spans the redshift range (0.1,1.8) has been suitably homogenized. 
Using the following merit function : 



we obtain thaix^-ed - 1-1 for 248 data points, and the best fit values are Hq - .91^^^^^y s - -1.49^qq4. In Fig (14) we compare 
the best fit curve with the observational data set, as in Paper 1. 
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Fig. 14. Updated Daly & Djorgovski database (Daly & Djorgovski 2005) fitted to our model. The solid curve is the best fit curve withx^^j 



1.1 



for 248 data points, and the best fit values are Hq = .91* 
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Fig. 15. We compare the sensitivity of the ^ relation to the values of the parameters in our model. Actually we plot the relative variation in ^ 
with respect to a variation of s from -1.4 to -1.3, and with respect to a variation of Ho from 1 to .9, the other parameters being fixed. The red 
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5.1 .2. H(z) and the relative galaxy ages 

In this section we discuss a possible observational determination of H{z) based on the method developed by 
(Jimenez & al. 2003)), that involves differential age measurements. We present some constraints that can be placed on the evo- 
lution of our quintessence model by such data. First, it is worth to stress some aspects connected with the sensitivity of the 
cosmology to the t{z) and ^ relations. Actually, it is well known that in scalar tensor theories of gravity, as well as in general 
relativity, the expansion history of the universe is determined by the function H(z). This implies that observational quantities, like 
the luminosity distance, the angular diameter distance, and the lookback time, all depend on H(z). It turns out that the most ap- 
propriate mathematical tool to study the sensitivity of the cosmological model to such observables is the functional derivative of 
the coiTesponding relations with respect to the cosmological parameters (see Saini, Padmanabhan & Bridle 2003 for a discussion 
about this point in relation to distance measurements). However, also from an empirical point of view, it is possible to show that 
the lookback time is much more sensitive to the cosmological model than other observables, like the luminosity distance, and the 
distance modulus. This circumstance encourages us to use, together with other more standard techniques discussed above, the 
age of cosmic clocks to test alternative cosmological scenarios. Apart from the advantage of providing an alternative instrument 
of investigation, the age-based methods use the stronger sensitivity to the cosmological parameters of the ^ relation, as shown 
in Figs. (15) and (16). Moreover, as we will discuss in the following, such a method reveals its full strength when applied to old 
objects at very high z. Actually it turns out that this kind of analysis could remove, or at least reduce, the degeneracy which we 
observe at lower redshifts, for example the one in the Hubble diagram for SNIa observations, which can be fitted by different 
cosmological models with similar statistical significance. Since the Hubble parameter can be related to the differential age of the 
universe as a function of redshift by the equation 



H(z) 



1 dz 
l+zdt' 



(77) 
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Fig. 16. We compare the sensitivity of the /j = in - M relation to the values of the parameters in our model in. Actually we plot 
the relative variation in fi with respect to a variation of s from -1.4 to -1.3, and with respect to a variation of Hq from 1 to 
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As wee see for the same variation in the parameters the sensitivity of the distance 



modulus is quite smaller than the respective variation in ^ 
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Fig. 17. The best fit curve of the measured values of H(z) corresponds to //q = 1.01 
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a determination of dz/dt directly measures H(z). Jimenez et al. ( Jimenez & al. 2003)) demonstrated the feasibility of this method 
by applying it to a z ~ sample of galaxies. With the availability of new galaxy surveys it becomes possible to determine H{z) 
at z > 0. Here we use the dz/dt data from (Simon & al. 2005)) to determine H(z) in the redshift range 0.1 < z < 1.8. To follow 
the procedure described in (Simon & al. 2005)), first we group together all galaxies that are within Az = 0.03 of each other This 
gives an estimate of the age of the universe at a given redshift. We then compute age differences only for those bins in redshift 
that are separated by more than Az = 0.1 but less than Az = 0.15. The first limit is imposed so that the age evolution between 
the two bins is larger than the error in the age determination. We note here that differential ages are less sensitive to systematic 
errors than absolute ages (see Jimenez & al. 2003) ). The observational value of H{z) is then directly computed by using Eq. (77), 
and after that the ages data have been scaled according to our choice of the unit of time (the unknown scaling factor has been 
provided by the;if^ procedure). To determine the best fit parameters, we define the following merit function : 



_ 1 



H{zr,s,Ho)-Hi 



(78) 



We obtain i\\?A.xled ~ ^-^^ ^ ^^^^ points, and the best fit values are Hq - 1.01+| 
the best fit curve with the observational data set. 



0.03' 



= -1.49 



+0.03 
-0.09- 



In Fig (17) we compare 



5.2. The Sunyaev-Zeldovich/X-ray method 

In this section we discuss how the parameters of our model can be constrained by the angular diameter distance as measured 
using the Sunyaev-Zeldovich effect (SZE) and the thermal bremsstrahlung (X-ray brightness data) for galaxy clusters. The dis- 
tance measurements using Sunyaev-Zeldovich effect and X-ray emission from the intracluster medium have to take into account 
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Fig. 18. Tile best fit curve to the H{z) data for our nmc model (dark blue line) and for the quintessence QCDM fitted to the new released 
WMAP- three years + SNLS data (WMAP New Three Year Results 2006) £1a = 0.72 ± 0.04, Sit = -O.OIO^JJ;™^, w = -l.Oe^SJi^^ (blue line). 
It is interesting to note that future high redshift data could disentangle the degeneration among different models, since at high z the predicted 
H(z) more sensitively depends on the values of the parameters. 



that these processes depend on different combinations of some of the parameters of the clusters (see Birkinshaw 1999 and refer- 
ences therein). The SZE is a result of the inverse Compton scattering of the CMB photons on hot electrons of the intercluster gas, 
which preserves the number of photons, but allows photons to gain energy and thus generates a decrement of the temperature in 
the Rayleigh- Jeans part of the black-body spectrum while an increment appears in the Wien region. We limit our analysis to the 
so called thermal or static SZE. The kinematic effect, present only in clusters with a nonzero peculiar velocity with respect to the 
Hubble flow along the line of sight, will be neglected since typically the thermal SZE is an order of magnitude larger than the 
kinematic one. As in Paper I, we introduce the so called Compton parameter, y, defined as the optical depth t - crj j rigdl times 
the energy gain per scattering: 



y '■ 



I 



n^crjdl. 



(79) 



where is the temperature of the electrons in the intracluster gas, m^ is the electron mass, is the number density of the 
electrons, and crj- is the Thompson cross section of electron scattering, and the integration is performed along the fine of sight. 
In the low frequency regime of the Rayleigh- Jeans approximation the shift of temperature is 



To 



(80) 



where Tq is the unperturbed CMB temperature. The next step to quantify the SZE decrement is to specify the model for the 
intracluster electron density and temperature distribution, which appropriately describes the observational properties of the gas. 
Following (Bonamente & al. 2005) we use a hydrostatic equilibrium double yS-model model. ActuaUy, at the center of clusters 
the density may be high enough that the radiative cooling time-scale is less than the cluster's age, leading to a reduction in 
temperature and an increase in central density. This can increase the central X-ray emissivity. At large radii, the density of the 
gas is sufliciently low that X-ray emission can be sustained for cosmological periods without significant cooUng. Therefore, cool 
core clusters effectively exhibit two components: a centrally concentrated gas peak and a broad, shallower distribution of the gas. 
This phenomenon motivated the modelling of the gas density with a function of the form: 
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(81) 



The quantity is the central density, / governs the fractional contributions of the narrow and broad components (0 < / < 1), 
and Vci are the two core radii that describe the shape of the inner and outer portions of the density distribution and fi determines 
the slope at large radii (the same is used for both the central and outer distribution in order to reduce the total number of 
degrees of freedom) ^ . This shape generalizes the single y6-model profile, introduced by Cavaliere and Fusco-Femiano (1976) 
and commonly used to fit X-ray surface brightness profiles, to a double y6-model of the density that has the freedom of following 
both the central spike in density and the more gentle outer distribution. A double j6-model of the surface brightness was first used 



^ It is worth to note that there are also different hydrodynamical models which predict universal gas density and gas temperature pro- 
files that agree with the observations (see for instance the one illustrated in ( Komatsu & Seljak 2001), derived from the universal dark 
matter density profile, assuming that the gas density traces the dark matter density in the outer parts of halos, or the one introduced in 
(Rasia, Tormen, Moscardini 2004)). 
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Fig. 19. Observational SZE data fitted to our model with the best fit values Hq = l+jj "^, s = -1.49![5°g, ^ = ^.70 ± 0.05. The empty boxes 
indicate distance measurements for a sample of 44 mentioned clusters (see Birkinshaw 1999, Reese & al. 2002),while the filled diamonds 
indicate the measurement of the angular diameter distances from Chandra X-ray imaging and Sunyaev-Zel'dovich effect mapping of 39 high 
redshift clusters (Bonamente & al. 2005). 



by Mohr et al. (1999) to fit X-ray data of galaxy clusters; the density model of Eq. (81) was further developed by La Roque 
(La Roque et al.(2006)). The X-ray surface brightness is related to the gas density as 

^--i^KlW/"^^-'' '''' 

where z is the cluster redshift, A^e is the X-ray cooling function, and it is a function of plasma temperature and energy in the 
rest frame of the cluster, including contributions from relativistic electron-ion thermal bremsstrahlung, electron-electron thermal 
bremsstrahlung, recombination, and two photon processes. The cluster angular diameter distance Da = dl/dO, where 6 is the 
line-of-sight angular size, can be inferred with a joint analysis of SZE, taking advantage of the different density dependence of 
the X-ray emission and SZE decrement: 

Sx'^ JnlAeedl^ J n]KeeDAde (83) 

ATcus^^nJJl^^n.,T.DAde. 

Actually it turns out that 

Da - (84) 

Note that Da is proportional to AT^^j^ and T^^^ (since A^e oc T^^), so the distance determination is strongly dependent on the 
accuracy of the SZE decrement and X-ray temperature measurements. 

Recently distances to 18 clusters with redshift ranging from z ~ 0.14 to z ~ 0.78 have been determined from a likeli- 
hood joint analysis of SZE and X-ray observations (see Table 7 in Reese & al. 2002). We perform our analysis using angular 
diameter distance measurements for a sample of 83 clusters, containing the 18 above mentioned clusters, other 24 known pre- 
viously (see Birkinshaw 1999), and a recently released sample with the measurement of the angular diameter distances from 
the Chandra X-ray imaging and Sunyaev-Zel'dovich effect mapping of 39 high redshift clusters of galaxies (0.14 < z < 0.89 
) (Bonamente & al. 2005). The unprecedented spatial resolution of Chandra combined with its simultaneous spectral resolution 
allows a more accurate determination of distances. Let us consider the merit function of the form : 
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(85) 



Fitting the data we obtain that jt'^^^ = 1.2 for 83 data points, and the best fit values are Ho - 1^()Q3, s = - 1 .49;^j)jjy and 
h - 0.70 + 0.05. In Fig. (19) we compare the best fit curve with the observational SZE data. 



5.3. Gamma-Ray Burst Hubble Diagram 

Gamma-ray bursts (GRBs) are bright explosions visible across most of the Universe, certainly out to redshifts of z ~ 7 and likely 
out to z ~ 10. Recent studies have pointed out that GRBs may be used as standard cosmological candles (Ghirlanda & al. 2004, 
Friedman & Bloom 2005): actually it turns out that the energy released during bursts spans nearly three orders of magnitude. 
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Fig. 20. Observational Hubble diagram for the recent SNIa sample compiled by (Riess et al. 2006) (empty lozenges), and the GRBs data 
(Bloom, Frail & Kulkarni 2003) (empty boxes) fitted to our model. The solid curve is the best fit curve with, Hq = 0.98^oq3, s = -1.431^ jj^. 
The red dashed line corresponds to the standard ACDM model with Q./^ = 0.71 



and the distribution of the opening angles of the emission, as deduced from the timing of the achromatic steepening of the 
afterglow emission, spans a similar wide range of values. However, when the apparently isotropic energy release and the conic 
opening of the emission are combined to infer the intrinsic, true energy release, the resulting distribution does not widen, as 
is expected for uncorrected data, but shrinks to a very well determined value (Frail & Kulkarni 2003), with a remarkably small 
(one-sided) scattering, corresponding to about a factor of 2 in total energy. Similar studies in the X-ray band have reproduced the 
same results. It is thus very tempting to study to what extent this property of GRBs makes them suitable cosmological standard 
candles. Schaefer (Schaefer 2003) proposed to use the two well known correlations of the GRBs luminosity (with variability, and 
with time delay), while others exploited the recently reported relationship between the beaming-corrected y-ray energy and the 
locally observed peak energy of GRBs (see for instance Dai & al. 2004). As for the possible variation of ambient density from 
burst to burst, which may widen the distribution of bursts energies. Frail & Kulkarni (Frail & Kulkarni 2003) remarked that this 
spread is already contained in their data sample, and yet the distribution of energy released is still very narrow. There are at least 
two reasons why GRBs are better than type la supernovae as cosmological candles. On the one hand, GRBs are easy to find 
and locate: even 1980s technology allowed BATSE to locate ~1 GRB per day, making the build-up of a 300-object database a 
one-year enterprise. The Swift satellite launched on 20 November 2004, detects GRBs at about the same rate as BATSE, but with 
a nearly perfect capacity for identifying their redshifts simultaneously with the afterglow observations Second, GRBs have 
been detected out to very high redshifts: even the current sample contains several events with z > 3, with one (GRB 000131) at 
z = 4.5 and another at z = 6.3. This should be contrasted with the difficulty of locating SN at z > 1, and absence of any SN with 
z > 2. On the other hand, the distribution of luminosities of SNIa is narrower than the distribution of energy released by GRBs, 
corresponding to a magnitude dispersion crM(SN) — 0.18 rather than o-m(grb) = 0.75. Thus GRBs may provide a complementary 
standard candle, out to distances which cannot be probed by SNIa, their major limitation being the larger intrinsic scatter of 
the energy released, as compared to the small scatter in peak luminosities of SNIa. There currently exists enough information 
to calibrate luminosity distances and independent redshifts for nine bursts (Schaefer 2003). These bursts were all detected by 
BATSE with redshifts measured from optical spectra of either the afterglow or the host galaxy. The highly unusual GRB980425 
(associated with supernova SN1998bw) is not included because it is likely to be qualitatively different from the classical GRBs. 
Bursts with red shifts that were not recorded by BATSE cannot yet have their observed parameters converted to energies and 
fluxes that are comparable with BATSE data. For the present analysis we shall use a sample of GRBs that had their redshifts 
estimated (Bloom, Frail & Kulkarni 2003), as represented by empty boxes in Fig. (20), with the distance modulus yu, given by 
Eq. (72). 

To this aim, the only difference with respect to the SNIa is that we slightly modify the correction term of Eq. (73), and we 
take 



^Geff(t) 

^G,„ = 2.5y ■ (86) 

'» (InlO)Ge// 

We expect that y is of order unity, so that the G-correction would be roughly half a magnitude. We obtain = 1.1 for 24 data 
points, and the best fit value is//o = 0.98^qq3, s - -1.43^qq4, which are compatible with the SNIa results. We also find that 
J which appears in Eq. (86) is equal to 1.5. In Fig. (20) we compare the best fit curve with both the GRBs and the SNIa Gold 
Sample. 



http://swift.gsfc.nasa.gov 
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5.4. The gas fraction in clusters 

Measurements of the gas mass fraction in galaxy clusters have been proposed as a test of cosmological models (Allen & al. 2002). 
Both theoretical arguments and numerical simulations predict that the baryonic mass fraction in the largest relaxed galaxy clus- 
ters should not depend on the redshift, and should provide an estimate of the cosmological baryonic density parameter Q.b 
(Eke & al. 1998). The baryonic content in galaxy clusters is dominated by the hot X - ray emitting intra-cluster gas so that what 
is actually measured is the gas mass fraction fgas and it is this quantity that should not depend on the redshift. Moreover, it is 
expected that the baryonic mass fraction in clusters is equal to the universal ratio Q.i,/Slm so that fgas should indeed be given 
by b-i£lb/Clm), where the multiplicative factor b is motivated by simulations that suggest that the gas fraction is lower than the 
universal ratio. Following the procedure described in (Allen & al. 2002, Allen & al. 2004), and already used in Paper 1 we adopt 
the standard CDM model (i.e., a flat universe with Q„ = 1 and h = 0.5) as a reference cosmology in making the measurements 
so that the theoretical expectation for the apparent variation of fgas with the redshift is: 



fgasiz) = 



(1 + 0.19 ^/h)^l„ 



(87) 



where we substitute the appropriate expression of Clm for our model, and D^'^^^ and D'^"'' are the angular diameter distance 
for the SCDM and our model respectively. Allen & al. (Allen & al. 2002) have extensively analyzed the set of simulations 
in (Eke & al. 1998) to get b = 0.824±0.089, so in our analysis below, we set b = 0.824. ActuaUy, we have checked that, for 
values in the 2cr range quoted above, the main results do not depend on b. Moreover we have defined the following merit 
function : 



I gas 

where 

i=l 



Qb/z^ -0.0214f (h-0.12f (b-0.824f 



0.0020 / ^\ 0.08 / ^\ 0.089 



fgas(zi,auHo)-f"^:(zdf 



(89) 



Here fgas(Zi) is the measured gas fraction in galaxy clusters at redshift z, with an error o-g, and the sum is over the Ngas clusters 
considered. Let us note that recently Allen & al. (Allen & al. 2004) have released a catalog of 26 large relaxed clusters with a 
precise measurement of both the gas mass fraction fgas and the redshift z. We use these data to perform our UkeUhood analysis, 
we getx^ = 1.16 for 26 data points, and s = -1.39+[;{j|, Hq = 0.98 ± 0.04, h = 0.65 ± 0.05, and = -0.95 ± 0.07. A brief 
comparison of our results with similar recent results of Lima et al. (Lima & al. 2003), where the equation of state characterizing 
the dark energy component is constrained by using galaxy clusters x-ray data, can still be done. In their analysis, however, they 
consider quintessence models in the standard gravity theories, with a non evolving equation of state, but they allow the so- 
called phantom dark energy with w < -1, which violates the null energy condition. As the best fit value of w to the data of 
(Allen & al. 2002) they obtain w = -1.29;|;Q7g2- In order to directly compare this result with our analysis we first fit the model 
considered in (Lima & al. 2003) to the updated and wider data set of (Allen & al. 2004), used in our analysis. To this aim we also 
refer to the model function fgasiz), and the merit function defined in the Eqs. (87) and (88) respectively. We gctx^ - 1-175 
for 26 data points, and Q„ = 0.23;|;gQ3, h = 0.76^^'^, and w = -1.11 ± 0.35, so w < -1, what corresponds to a phantom energy. 
Let us note that our model, instead, gives = -0.87 ± 0.05, what does not violate the null energy condition. In Fig. (21) we 
compare the best fit curves for our and the Lima & al. model with the observational data. 



6. Growth of density perturbations 

In this section we consider the evolution of scalar density perturbations in the longitudinal gauge ds^ = -(I + 2^)dt^ + a^(l - 
2(^)dx^. Contrary to what happens in the framework of the minimally coupled theory, where we have to deal with a fully 
relativistic component, which becomes homogeneous on scales smaller than the horizon - so that standard quintessence cannot 
cluster on such scales - in the non minimally coupled quintessence theories it is possible to separate a pure gravitational term 
both in the stress-energy tensor T^y, and in the energy density p^, so the situation changes, and it is necessary to consider also 
fluctuations of the scalar field. However, it turns out (Boisseau & al. 2000 , Riazuelo & Uzan 2002) that the equation for dust hke 
matter density perturbations inside the horizon can be written as follows: 

6„ + 2H6„ - iccav P,n S„, - , (90) 



where Gcm is the effective gravitational constant defined by 



1 



1 + I6^(s) 



2^{sW 1 + 12^(5) 



(91) 
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Fig. 21. In tlie diagram we plot tlie best fit curve to the /ja, data for our nmc model (red thick line) and for the quintessence model (black thick 
line) considered in (Lima & al. 2003). It is interesting to note that, as pointed out also for the model described in (Demianski & al. 2006), even 
if the statistical significance of the best fit procedure for these two models is comparable, the best fit relative to our nmc model seems to be 
dominated by smaller redshift data, while the one relative to the Lima & al. model by higher redshift data. 



The equation (90) describes, in the non minimally coupled models, evolution of the CDM density contrast, 6„, = 6pmlpm, for 
perturbations inside the horizon. In our model the Eq.(90) is rather complicated and takes the form 



6m + 2 



fA-24((,(2^ + 9) + 6)B(^)f^ +(2^(^ + 3) + 3) ((;^)"' - B(,)))((^)"' + (f^ - l)z?(^)) 



A{s) 



AUs + l)(s + 2)ffe^ ((^)"' + B(s) {-^ 



s+r)(s+2)Vt~ 



2.S+3 



2C(s)\2s + 3)2 



1 + 16^(5) 

1 + I2^(s) 



■ 



(92) 



Eq. (92) does not admit exact solutions, and can be solved only numerically. However, since with our choice of normalization 
the whole history of the Universe is confined to the range f € [0, 1] and therefore to study the behavior of solutions for f ^ 
we can always expand the functions in Eq. (92) in series around f = 0, in order to get approximate solutions. Actually we obtain 
an integrable Fuchsian differential equation, which is a hypergeometric equation. We then use such a solution to set the initial 
conditions at f = to numerically integrate Eq. (92) in the whole range [0, 1]. We use the growing mode 6+ and define the growth 
index / as 

/.^, (93) 
a In a 

where a is the scale factor. Once we know how the growth index / evolves with redshift and how it depends on our model 

parameters, we can use the available observational data to estimate the values of such parameters, and the present value of Q^. 

The 2dFGRS team has recently collected positions and redshifts of about 220000 galaxies and presented a detailed analysis of the 

/ 

two-point correlation function. They measured the redshift distortion parameter yS = — , where b is the bias parameter describing 

b 

the difference in the distribution of galaxies and mass, and obtained that yS|7^o.i5 = 0.49 + 0.09 and b - 1.04 + 0.11. From the 
observationally determined p and b it is now straightforward to get the value of the growth index at z = 0.15 corresponding 
to the eff'ective depth of the survey. Verde & al. (Verde & al. 2001) used the bispectrum of 2dFGRS galaxies, and Lahav & al. 
(Lahav & al. 2002) combined the 2dFGRS data with CMB data, and they obtained 



bverde = 1.04+0.11, 

bi,H,,. - 1.19 + 0.09. 

Using these two values for b we calculated the value of the growth index / at z = 0.15, we get respectively 

/i = 0.51+0.1 , 
fo = 0.58 + 0.11 . 



(94) 
(95) 



(96) 
(97) 



To evaluate the growth index atz = 0.15 we first have to invert the z - f relation and find f(0.15): actually the z - f relation is rather 
involved and cannot be analytically inverted, so we perform this inversion numerically. Finally, we get s = -1.4 + 0.1, //q = 




Fig. 22. The growth index / in different cosmological models. The thick dashed red line corresponds to our non minimally coupled model. The 
blue thin dashed curve corresponds to the standard A-CDM model with n„, = 0.25, and the black solid line corresponds to another quintessence 
model with an exponential potential (described in Demianski & al. 2005). 

0.94;^g g2' ^0 = 0.5 + 0.06 which corresponds to ^A,ff = 0.65 + 0.08. In Fig. (22) we show how the growth index is changing 
with redshift in our non minimally coupled model as compared with the standard A-CDM model, with Q.,„ = 0.25, and a 
quintessence model namely the minimally coupled exponential model described in (Demianski & al. 2005). We note that at low 
redshift theoretical predictions of these different models are not distinguishable, independent measurements from large redshift 
surveys at different depths could disentangle this degeneracy. 

7. Summary and conclusions 

In this paper we have extended the analysis that we performed in Paper I (where we have analyzed a special extended quintessence 
model, based on one of the most commonly used quintessence potentials V{(p) - Aij/', corresponding to the coupling F{^) - 
(3/32)0^ ), considering a new and wider class of theories for which exact solutions of the Einstein equations are known, and 
discussed how in such models it is possible to treat the fine tuning problem in an alternative way. We have shown that in the 
family of such models selected by requiring that their corresponding point like Lagrangian admits a Noether symmetry an epoch 
of accelerated expansion appears in a natural way. In the non minimally coupled scalar tensor theory of gravity it is possible 
to perform an appropriate conformal transformation and to move from the Jordan picture to the standard Einstein one but then 
matter becomes coupled to the scalar field. We have explored both descriptions and also considered the neutrino mass varying 
model as a possible example of non minimally coupled scalar tensor theory. 

It turns out that the imposed requirement of existence of a Noether symmetry is quite restrictive and we obtained a family of 
models that is fully specified by 3 parameters: a parameter s that determines the strength of the non minimal coupling and the 
potential of the scalar field, Hq the Hubble constant, and a parameter Vq that determines the scale of the potential. To determine 
the values of these parameters we compared predictions of our model with several independent observational data. The results of 
this parameter determination procedure are presented in Table 2. We see that with our average value of s = - 1 .46 the scale factor, 
for small f, is changing as a ~ fi-^^ ~ , and for large f, as a ~ f^-^^ ~ /^^"^ while ~ f0.0032 ^j^^ ^ ^ ^0.083 -^^ corresponding 
asymptotic regimes. The potential V is decaying to zero for large f, after reaching a maximum value, (see Fig. (23)). Similarly, 
the effective gravitational coupling Geff - ~jf '^^ decreasing for large f, until it becomes zero for f — » oo ( we have a sort 
of asymptotic freedom at f — > oo). It turns out that in our model the observational constraints on the variation of the effective 
gravitational constant are respected. Actually a new analysis of the Big Bang Nucleosynthesis (Copi et al. 2004) restricts the 
variations of G to 

AGbb^/Go = 0.01!°:2« (98) 

at 68% confidence level, where AGbbn/Gq = 2^^^, and Go = Gi^ewton- The combined analysis of the new He and WMAP 
data implies that 

-0.10 <AGbba,/Go< 0.13. (99) 

A recent analysis of the secular variation of the period of nonradial pulsations of the white dwarf G117-B15A 
shows (Benvenuto ef fl/. 2004) that < G/G < 4.0 x 10 " yr"' at 2cr, which is of the same order of magnitude as previous 
independent bounds (see also Biesiada & Malec 2004). With our unit of time this becomes < G/G < 4.0 t x 10"^, where 
T = 13.73+" }^ as given by the WMAP team (Spergel & al. 2006). 
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Fig. 23. The potential V as a function of the redshift z- 

Initially, for small f, the matter energy density is larger than the energy density of the scalar field. 

In Table (7) we present results of our analysis, they show that predictions of our model are fully compatible with the recent 
observational data. Comparing results of this paper with our previous analysis of minimally coupled scalar field models (see 
Paper I) we conclude that the present day observational data connected with the post recombination evolution of the universe can 
be fitted by several different models of quintessence. More data on high redshift supernovae of type la and GRBs are needed as 
well as more information on the early phase of structure formation in order to place stronger restrictions on the allowed type of 
dark energy. 

Table 2. The basic cosmological parameters derived from our model are compared with observational data. 



dataset 


s 


Ho 






high redshift SNIa 


-1 46+''"-* 

^•^"-0.03 


"•'^ -0.05 


69+" "^ 


-1.01 ±0.02 


dimensionless coordinates test 


_i 49+0.02 


n 97+0.04 

"•^'-0.03 


0.73 ± 0.08 


-0.98 ± 0.03 


relative galaxy ages 


_| 4Q+0.03 


1 01+'"" 

^•"^^-0.03 


0.80 ± 0.07 


-0.99 ± 0.04 


sze data 


_i 49+0.03 


1 +0.01 
^-0.03 


0.70 ± 0.05 


-1 ±0.04 


GRBs dataset 


_i 4Q+0.02 
^•^-'-0.04 


"•^°-0.03 


0.76 ± 0.06 


-0.9 ± 0.04 


fraction of gas in clusters 


1 19+0.04 
^•-^^-0.01 


0.98 ± 0.04 


77+""'^ 

' ' -0.03 


-0.87 ± 0.05 


galaxies peculiar velocity 


-1.4±0.1 


"•=^^-0.02 


0.65 ± 0.08 


-0.85 ± 0.04 


average 


-1.46 ±0.09 


0.97 ± 0.02 


0.74 ± 0.03 


-1.01 ±0.02 
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